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QiQ i The >100 GeV y-ray source, HESS J1713— 381, apparently associated with the shell-type supernova remnant (SNR) CTB 37B, was discovered 
f*"^ , using H.E.S.S. in 2006. X-ray follow-up observations with Chandra were performed in 2007 with the aim of identifying a synchrotron counterpart 
to the TeV source and/or thermal emission from the SNR shell. These new Chandra data, together with additional TeV data, allow us to investigate 
the nature of this object in much greater detail than was previously possible. The new X-ray data reveal thermal emission from a ~4' region in close 
proximity to the radio shell of CTB 37B. The temperature of this emission implies an age for the remnant of ~5000 years and an ambient gas density 
of ~0.5 cirT 3 . Both these estimates are considerably uncertain due to the asymmetry of the SNR and possible modifications of the kinematics due 
■ to efficient cosmic ray (CR) acceleration. A bright («7 xlO~ 13 ergcrrT 2 s~') and unresolved (<1") source (CXOU J171405. 7-381031) with a soft 
(r « 3.3) non-thermal spectrum is also detected in coincidence with the radio shell. Absorption indicates a column density consistent with the 
thermal emission from the shell suggesting a genuine association rather than a chance alignment. The observed TeV morphology is consistent with 
an origin in the complete shell of CTB 37B. The lack of diffuse non-thermal X-ray emission suggests an origin of the y-ray emission via the decay 
of neutral pions produced in interactions of protons and nuclei, rather than inverse Compton (IC) emission from relativistic electrons. 
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1. Introduction 

Approximately 50 sources of very-high-energy (VHE; 
>100 GeV) y-rays have now been detected along th e plane of 
our Galaxy dAharonian et all l2006bt iHintonl l2007l) . Some of 
these sources are (at least partially) coincident with supernova 
remnants (SNRs), consistent with the paradigm of cosmic-ray 
acceleration at SNR shocks. Resolved VHE y-ray emission 
from the shells of two nearby [d~l kpc) SNRs has been 
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detected from: RX J0852.0-4622 dAharonian et all l2005h and 
RXJ1713.7-3946 dAharonian et all l2007l) . However, in more 
distant objects, emission from the shell is often not resolvable 
with current instruments and there is ambiguity in the origin of 
the VHE emissi on. In several cases pulsar wind nebulae (PWNe; 
for a review see iGaensler & Sland d2006l) ). powered by neutron 
stars left behind in the supernova responsible for the remnants, 
are viable alternative acceleration sites (for example in the cases 
of HESS J1813-178 dFunk et alj |2007ah. and HESS J1804-216 
dAharonian et al.ll2006bl) ). 

The sensitivity to detect much more distant y-ray SNRs 
(c/»lkpc), and hence learn more about the population of such 
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objects, has been achieved only in a few regions of our galaxy 
where deep observations have taken place with the most sensi- 
tive available TeV instruments. One such region is that within 
a few degrees of RXJ1713. 7-3946, where a >50hour obser- 
vation with the 4-telescope H.E.S.S. system results in a point- 
source sensitivity of 1.5 x 1CT 13 photons cirT 2 s _1 above 1 TeV 
(or ~0.7% of the flux of the Crab Nebula). Within this region 
(1 -348°, b -0°), a y-ray source, HESS J1713-381 was de- 
tected in coincidence wi th the CTB 37 region of radio emission 
dAharonian et alj l2006bl) . This region consists of three super- 
nova remnants, CTB 37A (which consists of two separate rem- 
nants overlapping in projection: G 348.5+0.1 and G 348.5+0.0) 
and CTB37B (G 348.7+0.3) dKassim et alj [1991) . The first 
deep TeV measurements resulted in a significant excess close 
to CTB 37B and a less significant excess close to CTB 37A 
dAharonian et alj [2006bl) . Including the more recent H.E.S.S. 
data discussed here, the latter source (HESS J1714-385) is de- 
tected with a high l evel of confidence and is discussed in detail 
in a separate paper ( Ah aronian et al.ll2008l) . 

At radi o wavelengths CTB 37A a nd B have very similar 
properties dClark & Stephenson1ll975b : similar s urface bright- 
nesse s (~0.2 Jy arcsec~ 2 ), spectral indices (a ~0.3 dKassim et alj 
1 19911) ) and sizes (~5'). CTB37A is thought to be interact- 
ing with surro unding molecular clouds ( with densities of 100- 
1000 cirT 3 ) dRevnoso & Man gum 2000) and appears to have 
signi ficant OH mase r emission, indicative of dense shock heated 
gas dFrail et alJll996l) . No maser emission has been confirmed 
from the direction of CTB 37B. 

The distance to CTB 37B has been estimated by a number of 
authors. Casw elTet al.1 d 19751) used HI absorption measurements 
to determine the kinematic distance to C TB 37B (using the rota- 
tion curve of ISchmidt & Blaauwl d 1965b ) to be 10.2+3.5 kpc. A 
similar result was found for the distance to CTB 37A leading to 
the claim that CTB 37A and B are close not only in projection 
but also along the line of sight. Reinterpreti ng this result using a 
more recent galactic rotation curve model dBrand & BlitzHl993l) 
yields a distance ra nge of 5-9 kpc (we note that the value adopted 
in the catalog ue of lGreenl d2006l) is 8 kpc). 

Th e age of CTB 37B is also uncertain. IClark & Stephenson] 
d 1975b give age estimates for both CTB 37A and B of 
~1500 years and claim that there is a possibility that either 
source could be the remnant of the supernova of AD 393. This 
claim is further investigated by Downe| d!984i) who concluded 
that although such an association is possible, the remnants may 
in fact be considerably older. 

Both to improve our understanding of the SNR CTB 37B and 
to investigate the nature of HESS J1713-381, sensitive X-ray 
measurements of this region were considered highly desirable. 
We observed the field of CTB 37B with Chandra in early 2007. 
These new X-ray data, together with additional H.E.S.S. data, 
are discussed here, with the aim of elucidating the relationship 
between these two objects. 

2. H.E.S.S. observations and analysis 

The region toward CTB 37 was observed between April and 
October 2004 using the High Energy Stereoscopic System 
(H.E.S.S.) during a survey of the Galactic Plane, and in pointed 
observations toward the SNR RX J 1 7 1 3 .7-3946, yield ing 37 
hours of on-source live-time (AharonianetaD |2006bj). VHE 
gamma-ray emission was detected at a 6.3<x significance level 
from a region of rms size 0.06°+ 0.04°, spatially coincident 
with the SNR complex, at a position RA = 17 h 13 m 57.6 s , 
Dec = 38°12'0" (epoch J2000) and was announced as a new 



source HESS J 17 1 3— 38 1 . Initial spectral analysis of this source 
led to a photon index of 2.27 + 0.48 and a flux above 200 GeV 
of (4.2 + 1.5) Xl0 _12 cm _2 s _1 . The region has since been re- 
observed with H.E.S.S. and the data set used for analysis here 
includes all telescope pointings within 1.8° from the best fit po- 
sition of the source, leading to a total on-source live-time of 55.4 
hours. Correcting this to the equivalent on-axis exposure yields 
40.8 hours. The observations were performed over a large range 
of zenith angles (14°-60°) with an average value of 35°. 

Analysis was carried out using stan dard H.E.S.S. procedures 
(for details see Aharonia n et al.1 d2006al) ). An independent model 
analysis dde Nauroisll2006f) was performed and yielded consis- 
tent results. Hard cuts (image size cut 200 photoelectrons (p.e.)) 
provide a higher signal to noise ratio and a better PSF for 2- 
D images at the expense of an increased energy threshold and 
were applied for morphological analysis. Standard cuts (80 p.e.) 
provide a wider energy range and were thus applied for spectral 
extraction. The average energy threshold of the data set was 370 
GeV for hard cuts and 250 GeV for standard cuts. 

Two different background estimation models were used in 
the an alysis of the data (described in detail in iBerge et aLI 
(2007)). For spectral analysis, a reflected-region background 
model was applied using only off-regions with angular displace- 
ment less than 1.8° from the best fit source position in an attempt 
to reduce systematic errors arising from the strong exposure gra- 
dient across the field of view (due to most pointed observations 
being directed at RX J1713. 7-3946). For morphological anal- 
ysis and 2-D image generation, used for determining the best 
fit position of the source, a ring background model was used 
with ring radius 0.5°, excluding all known y-ray sources such as 
HESS J1714-385 from the background region. 

The integration region used for estimating the flux and dif- 
ferential energy spectrum of the VHE emission was a circle of 
radius 0.187° centered on the best fit position. This size was 
chosen so as to maximise the encompassed excess whilst min- 
imising contamination from HESS J1714-385. From the source- 
fitting it was estimated that this region contained 93% of the total 
emission fromHESS J1713-381 and included an estimated 11% 
contamination from HESS J1714-385. 

Figure Q] shows a smoothed, background subtracted and ac- 
ceptance corrected image of the region surrounding CTB 37. 
TeV emission is observed from both HESS J1713— 381 (the 
centroid of which is marked with a filled black triangle) and 
HESS J1714-385. Radio contours, t aken from the Mo longlo 
Galactic Plane Survey at 843 MHz dGreen et alj|1999h . show 
emission only from the Eastern side of the SNR CTB 37B as 
well as a faint plateau of emission extending to the South-East 
of the remnant. A circle of radius 5.1' has been drawn onto 
the image in Figure Q] to represent the size of the SNR based 
on the partial radio shell. The centroid of the TeV emission 
from HESS J1713-381 is located at RA = 17 h 13 m 54 s + 4 s stat , 
Dec = -38°11'58" + 45" Jto/ , close to the centre of this pro- 
posed shell. The systematic error in the source location is esti- 
mated to be 20" in both coordinates. The source shows an ex- 
cess of 292 + 36 counts with a st a tistica l significance (using 
the likelihood method of lLi&Mal d 1983b ) of 8.6cr. A some- 
what higher significan ce is obtained using the model analysis 
dAharonian et alj|2008h . 

A power law with photon index T = 2.65 ± 0.19 JMf + 0.20 JVJ 
and a differential normalization at 1 TeV of 
(6.5 ± l.l stat + l^JxlO^cnrV'TeV" 1 provides a 
satisfactory description of the data (^ 2 /dof = 3.8/3) as shown 
in Figure [2] The integral flux above 200 GeV was found to be 
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Fig. 1. Acceptance corrected y-ray excess map of the SNR com- 
plex CTB 37 after Gaussian smoothing with a radius cr sm of 2.9'. 
The black contours represent 0.08, 0.16, 0.32, 0.64, 1.28 and 
2.56 Jy/beam radio emiss ion from the Molo nglo Galactic Plane 
Survey (beam size ~1') dGreen et al.| [T999). The field of view 
of the Chandra observation is shown by a white box surround- 
ing CTB 37B. The estimated si ze of the supe rnova remnant, cen- 
tered on the position stated by iGreerJ d2006f) . is illustrated by a 
dashed white circle. The simulated PSF is shown in the box in 
the bottom right corner and is smoothed in the same way as the 
excess map. The colour scale of the map is in units of counts per 
2ncr\ m = counts per 0.0145 degrees 2 . 
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Fig. 2. Differential energy spectrum of HESS J1713-381. Solid 
symbols are derived from the full data set and are fitted with a 
power law with photon index F = 2.65 ± 0. 19 (solid line). A 90% 
confidence upper limit is given at energies above 10 TeV. The 
empty squares and the dashe d line represent the data acquired 
during previous observations (Aharoni an et al ] l2006bh . 
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Fig. 3. X-ray image (1-5 keV) of the Chandra field of view 
showing the positions of the 18 detected point sources on the 
four ACIS CCD chips. The image is a count map after Gaussian 
smoothing with <x sm = 5". Source #1 is located very close to 
one of the CCD chip edges. The position of the centroid of the 
TeV emission from HESS J 17 1 3-38 1 is marked by a white tri- 
angle. The white box shows the region from which the emission 
profile in Figure [5] is taken and the dashed white ellipse shows 
the region from which the diffuse X-ray spectrum was extracted. 
The colour scale shows counts per 2n(j\ m and has been trun- 
cated to better show the fainter point sources in contrast to the 
background. The maximum prior to truncation was 1214 at the 
position of source #1. This truncation forces the brighter sources 
to appear larger than their true size. As in Figure Q] the dotted 
circle illustrates the estimated size of the SNR CTB 37B 



(5.6 ± \.5 sm ±l.l jyj )xl0 cm s and the energy flux in the 
range 0.5-5 TeV is ^2xl0 _12 erg cm _2 s _1 . 

A joint morphological fit of HESS J1714-385 and 
HESS J1713— 381 was made with a model consisting of 
two Gaussian emission profiles convolved with the PSF of 
the instrument. The intrinsic rms width of HESS J1713-381 
extracted from this fit, 2.6' + 0.8' , is compatible with the size 
found during the H.E.S.S. survey dAharonian et al 1 l2006bl) . An 
unresolved thin shell of radius r has an equivalent Gaussian rms 
of <x ^0.6r. For a filled sphere <x ^0.45r. The measured width 
is consistent with a shell of radius ~4-6', compatible with the 
estimated radio shell size. 



3. X-ray observations and analysis 

The CTB 37B region was observed with the Advanced 
CCD Imaging Spectrometer (ACIS) of the Chandra X-ray 
Observatory on the 2 nd February 2007 for 26 ks (obser- 
vation ID 6692). The raw data were analysed using the 
Chandra Interactive Analysis of Observations (CIAO ver- 
sion 3.4, CALDB version 3.4. 1). The data set was unaffected by 
soft proton flares enabling the full observation time to be utilised 
for analysis purposes. 

Figure [3] shows a count map of the region after Gaussian 
smoothing with a 5" radius to match the off-axis PSF. This im- 
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Fig. 4. Exposure corrected count map (1-5 keV) after subtraction 
of emission from the 18 point sources detected in the field of 
view. The image has been Gaussian smoothed with <x sm = 50" 
to highlight the diffuse emission. The colour scale is in units of 
counts per 2ncr\ m . The black circle shows the position of source 
#1. White square symbols show the location of the other point 
sources which have been removed from the image. The white 
contours represent the same radio contours as in Figure Q] The 
position of the centroid of HESS J1713— 381 is shown as a white 
triangle. 



age illustrates the positions of 18 point-like sources, described 
in Table Q] detected at a >5<x level within the Chandra field of 
view using the CIAO wavdetect algorithm. The most significant 
of these sources, source #1, is located 3' North-East of the cen- 
troid of the TeV source HESS J1713-381. 

Figure |4] shows an exposure corrected, smoothed excess map 
after removal of emission from the point sources described in 
Tablefl] Extended diffuse emission is observed at a > 1 1 cr level to 
the East of HESS J 17 13-381 and appears to be contained within 
the region defined by the radio contours. The diffuse emission 
extends over all four ACIS CCD chips. This peak is spatially 
coincident with the radio shell but is not aligned with the peak 
of the radio emission and is instead offset to the West by ~1.6'. 
This shift is highly significant considering the 0.6" astrometric 
accuracy of Chandra. No evidence for significant temperature 
variations across the diffuse emission was found and spectral 
properties on each individual ACIS chip were consistent. Hence 
the region used to extract the reported diffuse spectrum extends 
over all four chips, and is indicated by a dashed white ellipse in 
Figures [3] and |4] 

Figure [5] shows a profile of the X-ray emission in the 
rectangular region shown in Figure [3] after the removal of 
the 18 point sources found within the Chandra field of view. 
This region is centered on the position of the X-ray source 
CXOUJ171405.7-381031 (source #1). The simulated PSF 
model for source #1 is shown in red. The rms width of the diffuse 
emission extracted from a Gaussian fit to the data (black curve, 
* 2 /dof = 86.6/57) is cr = 1.2' ±0.1'. 

The energy spectrum of source #1 was taken from an integra- 
tion region of radius ~ 1 .7" centered on the best fit position. This 



S 120 

^100 



80 



60 



40 



a. 20 



-20 



- ]> 










-0.08 -0.06 -0.04 -0.02 



0.02 0.04 0.06 0.08 

Relative Position (°) 



Fig. 5. Acceptance corrected emission profile taken from the 
rectangular region shown in Figure[3] after the removal of emis- 
sion from the 18 detected point sources. This region is centered 
on the best fit position of source #1. A Gaussian fit to the data is 
shown by the smooth curve. The simulated PSF model for (the 
subtracted) source #1 is shown in red. Source #1 is ~10 times 
brighter (~1200 counts) than the observed diffuse emission. The 
background is estimated from regions of equal size and shape 
displaced to the north and south of the diffuse emission. 



spectrum is well described by a single absorbed power law with 
N H = (4.1+°*) xl0 22 cm -2 and r _ 3 3+0.2 with x 2 /do{ _ 68 5/ 6 3 



(the probability of a blackbody model fit to this source was 
~0.2%). Both sources #2 and #3 are consistent with blackbody 
models with kT<0.5 keV. The other 15 point sources were too 
weak to allow a meaningful spectral analysis to be performed. 

Regions around the 18 detected X-ray point sources were 
excluded from the data set before extracting the spectrum of the 
extended diffuse emission. The background was estimated from 
a large region covering all 4 ACIS chips and excluding regions 
around the point sources and the diffuse emission. The resulting 
spectrum was best fit with an absorbed non-equilibrium ionisa- 
tion (NEI) model with r = (7^) xl0 1() s cirT 3 0t 2 /dof = 68.9/62) 
and is shown in the bottom panel of Figure [6] A marginally ac- 
ceptable fit is also possible with t set to the maximum allowed 
value of 5 x 10 13 s cm -3 , indicating that the source is approach- 
ing equilibrium. No evidence was found for diffuse emission 
above 5 keV. Parameters of fits with different spectral models 
are given in Table [2] The consistent column densities from the 
model fits suggest a source #1 and the diffuse emission lie at 
approximately the same distance. 

Timing analysis of the source was difficult due to its po- 
sition on the chip edge leading to an artificial modulation at 
10~ 3 Hz caused by instrumental dither effects, and the 3.2 s frame 
time. However, we find no evidence for variability on 3xl0~ 3 to 
0.2 Hz timescales. 



4. Discussion 

4.1. The X-ray emission of the SNR CTB 37B 

The X-ray measurements presented here allow a better as- 
sessment of the age and environment of CTB 37B than was 
previously possible. From the age and distance of the remnant an 
external density can be estimated. Unfortunately the distance to 
the remnant is rather uncertain. Here we adopt a 7 kpc distance 
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ID 


Name 


RA 


Dec 


Counts 


Significance 




CXOUJ 


17 h 


-38° 




cr 


1 


171405.7-381031 


14 m 5.758 s 


10'31.32" 


1171 


158.9 


2 


171354.1-381844 


13 m 54.129 s 


18'44.24" 


445 


38.6 


3 


171404.3-381237 


14 m 4.385 s 


12'37.51" 


64 


23.2 


4 


171421.3-381310 


14 ra 21.38 s 


13' 10. 19" 


55 


19.6 


5 


171438.4-381114 


14 ra 38.43 s 


11' 14.2" 


40 


11.4 


6 


171323.4-381641 


13 m 23.45 s 


16'41.8" 


40 


9.8 


7 


171424.4-380959 


14 m 24.49 s 


9'59.1" 


28 


9.1 


8 


171402.2-381115 


14 m 2.224 s 


11' 15.4" 


19 


8.4 


9 


171409.7-380652 


14 m 9.77 s 


6'52.7" 


23 


8.2 


10 


171404.0-381652 


14 m 4.09 s 


16'52.4" 


23 


7.6 


11 


171444.3-380910 


14 m 44.35 s 


9' 10.5" 


19 


6.9 


12 


171346.0-381705 


13 m 46.08 s 


17'5.8" 


37 


6.5 


13 


m A 1 1 ^ Tom c 
1 / 1413.6-381215 


14 13.66 


12 15. / 


14 


6.4 


14 


171436.8-381331 


14 m 36.88 s 


13'31.8" 


31 


6.3 


15 


171341.0-381558 


13 m 4 1.07 s 


15'58.7" 


19 


5.9 


16 


171339.6-381320 


13 m 39.61 s 


13'20.0" 


15 


5.4 


17 


171344.2-380822 


13 m 44.22 s 


8'22.0" 


13 


5.1 


18 


171428.2-380545 


14 m 28.27 s 


5'45.6" 


23 


5.1 



5 



Table 1. Properties of all point sources detected above a 5<x level in the Chandra field of view using the CIAO wavdetect algorithm. 
Statistical errors on the positions are ~3" for the weakest sources on each axis. The astrometric accuracy of Chandra is 0.6" 



Element 


Model 




Best Fit Parameters 




Flux (1-5 keV) 


X z /dof 




(with absorption) 


N H (10 22 cm" 2 ) 


T kT (keV) 


T(10'°scm- 3 ) 


(10~ 13 erg cm^s" 1 ) 




Source #1 


Power law 


4 l +u - 3 


-j "3+0.2" 
■'"'-0.1 




6.7 


68.5/63 


Source #1 


Blackbody 


2.1 ±0.2 


0.76 ± 0.03 




5.9 


100.2/63 


Source #1 


Two-temperature 


9 q+0.2 
z..7_ 03 


52 +0 06 




6.4 


68.0/61 


Diffuse emission 


Blackbody 
NEI 


-i q+0.3 
J - -0.2 


i.u _ Q3 

0.8 ±0.1 


■7+9 


2.9 


68.8/62 


Diffuse emission 


NEI 


t O+0.5 
3.8_ 04 


0.7 ±0.1 


5xl0 3 (fixed) 


2.7 


75.5/63 


Diffuse emission 


Power law 


4 7+0-6 


r q+0.6 
J - J -0.5 




2.8 


101.6/63 



Table 2. Best fit spectral parameters for the point source CXOU J171405. 7-381031 and the diffuse X-ray emission with various 
models. Errors are given at the lcr level. The fluxes listed are the absorbed fluxes 



based on the HI absorption measurements of ICaswell et al.l 
(119751) scaled to more recent galactic rotation measurements. 
Assuming that the thermal X-ray diffuse emission is indeed 
associated with the SNR and that electrons and ions are in 
thermal equilibrium, the temperature obtained from the NEI 
fit (kT = 0.8 ± 0.1 keV) can be used to estima te the SNR age 
using the Sedov-Taylor phase equations (e.g. iBorkowski et al.l 
1200 lb . The temperature of the thermal emission implies a 
shock velocity, v s ~800(r/0.8keV) I/2 Ou/0-62m H ) _1/2 kms- 1 
and the age of the remnant is estimated to be 
-4900 (d/7kpc)(yu/0.62m H ) 1/2 (r/0.8keV)- I/2 years (where 
jj is the mean mass per particle). At a distance of 7 kpc the 
estimated radius of 5' corresponds to an intrinsic size close to 
10 pc. A spherical Sedov expansion then implies an external 
density of ~0.8 (Ea/10 51 erg) cm -3 where is the kinetic 
energy of the SN explosion. An age of 1600 years, compatible 
with an origin in the supernova of AD 393, would imply 
n -0.1 (Ei/10 51 erg) cm 3 . 

There are three main complications which could affect these 
estimates. Firstly, the Sedov solution holds only approximately, 
as the region surrounding CTB 37B appears to be inhomoge- 
neous, leading to an asymmetric expansion. Secondly, the tem- 
perature equilibrium between electrons and ions as s umed above 
is not seen in many young SNRs, see iRakowskH d2005l) for a 
recent review. In the case of non-equilibration, the velocity es- 
timate must be increased by a factor V(l + Tj/T e )/2 (asum- 
ing a pure hydrogen gas), making the remnant younger by the 
same factor. Finally, efficient acceleration of cosmic rays at the 



shock would result in less heating of the thermal gas for a given 
shock velocity (because a fraction of the dissipated kinetic en- 
ergy goes into the non-thermal population). The more efficient 
the acceleration, the faster the shock has to be to achieve a given 
gas temperature, leading to an over estimate of the SNR age 
using the conventional Sedov-Taylor gas dynamic solution. In 
the approximation of a s teady-state, plane-parallel shock (see 
e.g. iDrurv & Voelld[l98ll) . it may be shown that if the fraction 
of the downstream enthalpy in accelerated particles is then 
the velocity estimate must be increased by a factor of approx- 
imately 1 / Vl - and the age estimate decreased accordingly. 
For example, the combined effect of non-equilibrium tempera- 
tures, Ti - 3T e , and very efficient CR acceleration, # =0.4, 
would change the estimated age of the remnant to ~2700 years. 

The thermal X-ray flux can be used as a cross-check 
for this scenario. Due to absorption the luminosity im- 
plied by the measured flux is only a small fraction of the 
radiated power. The de-absorbed bolometric flux, fi, i, is 
*2.6xl0- n erg cm 2 s -1 . At 0.8 ke V (~10 7 K) the cooling coeffi - 
cient, A is «2.7xl0~ 23 erg cm 3 s -1 dSutherland & Dopital[l993l) . 
The ion density in the emission region is connected to fb i via 
fboi = A V n 2 1 4n d 2 , assuming thermal equilibrium and equal 
number densities for electrons and ions. Taking the volume, V, 
to be that of an ellipsoid with principle axes equal to the FWHM 
of the X-ray emision along its major axis, a = 2.8' (as shown in 
Figure |5), the FWHM in the perpendicular direction, b — 1 .9', 
and the FWHM along the line of sight (estimated as / X a), 



6 



H.E.S.S. Collaboration: Chandra and H.E.S.S. observations of the SNR CTB 37B 



~ 0.02 



0.015 



0.01 



0.005 




0.04 



0.03 



0.02 



Diffuse Region 




6 7 
Energy (keV) 



Fig. 6. Top: Count rate versus energy for source #1 (black filled 
symbols) compared to the best fit model for an absorbed power 
law. Bottom: Count rate versus energy for the diffuse emission 
compared to the best fit model for an absorbed NEI model. Fit 
parameters are described in the text. 



we find n ~2 (fd/7 kpc) -1 / 2 cirT 3 . The non-observation of ther- 
mal X-ray emission in other parts of the shell implies n < 1 
cirT 3 in these regions. If this emission occurs in a shock com- 
pressed region the external density is likely a factor ~4-7 lower 
and hence broadly consistent with the previous estimate. A gen- 
eral picture therefore emerges of CTB 37 B as a few 1000 year 
old SNR expanding into a (somewhat inhomogeneous) medium 
with average density n ~0.5 cm 4 . This discussion assumes that 
the bulk of the thermal X-ray emission originates in the shocked 
surrounding medium; if a significant fraction were to originate 
in shocked SN ejecta, the above estimates would have to be sub- 
stantially revised. 



shell, it seems plausible that there is a physical association of 
these two objects. In this context, source #1 may represent the 
X-ray emission of an isolated neutron star. We find no clear 
counterpart for this object at optical/infrared wavelengths; the 
closest catalogued source is the star GSC 2.3 S8PL008678 

from the Guide Star Catalogue (Version 2.3.2) 

( Spac e Telescope Science Institute & Osservatorio Astronomico di Torino! 
120071). offset by ~4~ The- closest 8/im GLIMPSE 



(Benjamin et al. 2003) source apparent by visual inspection lies 
3" away. 



Whilst the luminosity of source #1 (~4 x 10 erg s ) 
lies in a range fairly typical for the X-ray nebulae of Vela- 
like PWN (for example PS RJ1823-13 and PSRJ181 1-1925, 
see iKargaltsev et al.l (120071) ). the lack of extended, non-thermal 
emission around source #1 argues strongly against the PWN hy- 
pothesis. In addition, an identification of source #1 with the neu- 
tron star from the SN explosion of CTB 37B would imply a large 
proper motion of this object. Assuming a birth place of the pul- 
sar at the geometric centre of the remnant (RA = 17 h 13 m 55 s , 
Dec = -38°11'0") dGreenl 120061) and an age t yields an im- 
plied 'kick' velocity of ~1000(c//7kpc)(f/5000 years)" 1 kms" 1 . 
This estimate lies towards the upper edge of th e range of veloc- 
ities observed for pulsars (10-~2000kms -1 ) dj ben & Tutukov 
119961: IChatterjee & Cordesll2002t lArzoumanian et al.1120021] 1 !. A 
pulsar with such a velocity is likely to be associated with a bow- 
shock morphology pulsar wind nebula, as seen for example in 
G 359.23-0.82 ("the Mouse") dGaensler et al.ll2004 where the 
pulsar has a velocity of a^OOkms -1 . No such feature is detected 
in this case. 

Source #2 is spatially coincident with the star 
USNO-A2.0 0450-26795004 and has a soft thermal spec- 
trum (kT ~ 0.2 keV) and relatively low N H of <6xl0 21 cirT 2 , 
consistent with emission from a nearby star. Source #3 also 
has a low temperature (~0.4 keV) black-body spectrum but the 
closest catalogued star is ~3" away. This source is spatially 
coincident with the diffuse X-ray emission, and has a consistent 
column density ((6±2) xl0 22 crrT 2 ) suggesting that these objects 
may be physically associated. The closest source to the centroid 
of VHE emission is source #8 but the 19 excess counts are not 
sufficient for spectral analysis. This source has no catalogued 
stellar counterpart. 



4.3. The nature of the y-ray source HESSJ1 713-381 

The position and size of HESS J 17 13—38 1 make an association 
with CTB 37B very likely. Ho wever, as in the ca se of other ob- 
jects such as HESS Jl 8 13- 178 dFunk et alj|2007al) . there is pos- 
sible ambiguity between an origin of the TeV emission in the 
SNR shell and a pulsar wind nebula interpretation. Here we dis- 
cuss these alternative explanations in turn. 



4.2. X-ray point sources 

None of the 18 detected X-ray point sources were found to 
be spatially coincident with the centroid of the observed VHE 
emission and there is no obvious candidate for a compact central 
object of the SNR CTB 37B. For the three brightest sources (#1, 
#2, and #3) there are sufficient statistics for spectral analysis 
to be attempted and all are apparently point-like in nature. 
Source #1 (CXOU J171405.7-381031) is by far the brightest 
and lies in coincidence with the radio shell of CTB 37B and 
close to the maximum of the diffuse X-ray emission of the 
shell. As it has a measured Nh consistent with that of the SNR 



4.3.1. TeV emission from a PWN 

There are now >10 associations of young pulsars and TeV y- 
ray sources. These associations suggest a typical fraction of 
~1% for the conversion o f rotational energy into TeV emis- 
sion dCarrigan et al.l l2007l) . Assuming that the explosion of 
CTB 37B left behind a pulsar, the TeV emission could be ex- 
plained as inverse Compton emission from "relic" electrons ac- 



1 We note that both the X-ray and radio emission suggest a denser 
medium on the Eastern side of the remnant implying that the explosion 
center may have been significantly further East than the current geomet- 
rical center, hence reducing the required proper motion of the pulsar 
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celerated soon after the pulsar birth. Assuming IC scattering 
only on the cosmic microwave background (CMB), the syn- 
chrotron cooling time of electrons responsible for TeV emis- 
sion is -8000 (B/10jUGr 2 (£ T /lTeV)- 1/2 years, long enough 
that electrons injected immediately after the pulsar birth would 
continue to contribute to the TeV emission. In contrast, 
~3 keV synchrotron emission is associated with higher en- 
ergy (~80(Z?/10//G) _1 / 2 TeV) electrons which cool somewhat 
faster: f sync * 1500(fi/10//Gr 3/2 (e sync /3keV)- 1/2 years, and 
would hence be confined to a region close to the present 
pulsar position. The protot ype of the class of TeV bright 
PWN, HESSJ1825-137 (lAharonian et alj l2006d) . is associ- 
ated with the ~2xl0 4 year old, high-spindown luminosity 
(E * 3xl0 36 ergs" 1 ) pulsar PSRB1823-13. HESS J1713-381 
has similar spectral shape (F = 2.65 + 0.19 versus F * 2.4) 
but is apparently somewhat smaller (diameter 10 pc cf 40 pc), 
younger (5 kyr cf 20 kyrs) and dimmer (0.5-5 TeV luminosity 
-l^xlO^ergs" 1 cf -8x10 s4 ergs" 1 ). 

However, even if CXOU J17 1405.7-38 103 1 is indeed a high 
kick velocity pulsar, the absence of extended non-thermal X-ray 
emission is very difficult to explain in this scenario. An explana- 
tion of HESS J 17 13-381 by inverse Compton emission from ul- 
tra relativistic electrons in the wake of the pulsar therefore seems 
rather implausible. 



4.3.2. TeV emission from the SNR shell 

Whilst there is a clear offset of the centroid of HESS J 17 13-381 
from the thermal X-ray and radio emission of the shell, the TeV 
centroid is consistent with the estimated geometrical centre of 
the remnant. In addition, the observed size of the TeV source 
is consistent with emission in a thin shell of radius ~4.5', con- 
sistent with the radius of the radio shell. If the observed radio 
and X-ray emission from the Eastern part of the shell is due 
to a local density and/or magnetic field enhancement, then it 
seems plausible that particle acceleration may also be taking 
place in the Western part. In a scenario where the TeV emis- 
sion originates from inverse Compton (IC) scattering of > TeV 
electrons, the emission is expected to directly trace the dis- 
tribution of such electrons (assuming uniform radiation fields 
on the scale of the SNR). However, X-ray synchrotron emis- 
sion is expected from the same population of electrons with 
a flux ~ F Te v(B /3/j.G) 2 /(U m d/ Ucmb), assuming that there is 
no cut-off in the electron energy spectrum between the range 
probed by the TeV IC emission and that probed by 2 - 10 keV 
synchrotron emission (see above). The lack of measured non- 
thermal X-ray emission in the shell implies an X-ray synchrotron 
flux F sync significantly lower than the observed thermal flux 
of ~3 x 10~ 13 erg ctrT 2 s _1 . Given the measured TeV flux of 
2x 10~ 12 erg crrT 2 s _1 (0.5-5 TeV), this IC scenario implies ei- 
ther an implausibly low B-field of ~ 1 fiG or a rather sharp cut-off 
in the electron spectrum around ~40 TeV. 

A more natural explanation may be an origin of the TeV 
emission via the decay of neutral pions produced in proton- 
proton interactions in the region of the SNR shell. Following 
Ahar onian et al.l (1 19941) the expected y-ray flux is F(>1 TeV) ~ 
lO- 12 0(£*/lO 51 erg) (d/7kpcT 2 (n/0.5 cm 4 ) cirT 2 s~' where 
is the CR acceleration efficiency. The measured flux above 1 TeV 
(3.9±0.7) xl0~ 13 cirT 2 s _1 implies efficient CR acceleration. 
The relatively soft measured TeV spectrum (F = 2.65 + 0.19) 
is somewhat surprising in this scenario but could be interpreted 
as the consequence of a cut-off in the proton spectrum at ener- 



gies somewhat lo wer than seen for exam ple in RX J 1 7 1 3 .7 -3946 
(£ max -100 TeV) (lAharonian et alj|2007l) . 



4.3.3. Association with the EGRET source 3EGJ1 71 4-3857 

HESS J1713-381 lies close to the 99% c onfidence contour of 
the >100 MeV source 3EG J1714-3857 dHartman et al.lll999l) 
and so this object could be considered as a potential counter- 
part to the observed TeV excess. The SNR RX 11713.7 -3946 has 
also b een suggested as a counterpart to this emission ( But t et all 
|2002|), as has the newly dis covered y-ray source HESS J1714- 
385 ( Aharonian et al.l 1200 8). Despite exhibiting much weaker 
TeV emission than RX J17 13.7-3946, HESS J1713-381 may be 
considered as a viable counterpart to 3EG J1714-3857 due to 
its softer spectral shape (F * 2.7). An extrapolation of the TeV 
spectrum down to GeV energies is marginally consistent with 
the measured flux of 3EG J1714-3857. However, such a spec- 
tral match is not unlikely for a chance pairing of H.E.S.S. and 
EGRET sources, given the respective instrum ental sensitivities 
and typical spectral shapes dFunk et aT1l2007b|). The angu lar res- 
olution and flux sensitivity of GLAST ( Thompsonll2004l) should 
allow this issue to be resolved in the near future. 

5. Summary 

The new X-ray and y-ray data presented here repre- 
sent a major improvement in our knowledge of the 
CTB 37B/HESS J1713— 381 system. The discovery of thermal 
X-ray emission from the shell allows the uncertainties on the age 
and ambient density of the remnant to be considerably reduced, 
with likely values of ~5000 years and ~0.5 cirT 3 emerging. The 
non-thermal point-source CXOU J171405.7-381031, which 
is detected embedded in the diffuse thermal emission, may 
be a neutron star associated with CTB 37B, but considerable 
difficulties exist with this interpretation. It seems likely that the 
TeV emission originates in the shell of CTB37B. Whilst an 
IC origin of the TeV emission cannot be excluded, the decay 
of neutral pions produced in proton-proton interactions in the 
(entire) shell of the supernova remnant seems to be the most 
natural explanation. 
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